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ABSTRACT 

We present an investigation of the X-ray absorption features predicted by hydrody- 
namic simulations of radio galaxies interacting with the intracluster medium (ICM) 
of their host galaxy clusters. We show how these absorption lines can be used as a 
new diagnostic for the radio-galaxy/ICM interactions. Such interactions have been 
observed in numerous systems by ROSAT, Chandra and XMM-Newton, and under- 
standing them has implications for AGN feedback and galaxy formation. Starting 
from the hydrodynamic simulations of Reynolds, Heinz & Begelman, we calculate the 
properties of the highly ionized iron and oxygen lines (seen in absorption against the 
central active galactic nucleus; AGN), predicting line shapes, equivalent widths, col- 
umn densities and velocity shifts. The main effect of the jet on the absorption lines is 
a reduction of the line strength from that of the quiescent ICM and the introduction 
of some velocity structure in the line profile. We investigate whether these features 
are detectable with current as well as future high-resolution X-ray spectrometers. We 
conclude that the Chandra transmission gratings have insufficient sensitivity to detect 
these features with high significance, and certainly would not allow a study of the 
dynamics of the interaction via absorption signatures. Constellation-X, on the other 
hand, will allow superb constraints to be derived. We can also use this analysis to as- 
sess the idea that radio-galaxy induced ICM outflows give rise to the resonant oxygen 
X-ray absorption lines that have been claimed as evidence for the warm-hot intergalac- 
tic medium (WHIM). We show that these detached, high- velocity oxygen absorption 
lines cannot result from a radio-galaxy/ICM interaction, thereby strengthening the 
WHIM interpretation. 

Key words: galaxies: jets — intergalactic medium — line: profiles — shock waves 
— X-rays:galaxies:clusters 



1 INTRODUCTION 

It has been known for over a decade that radio galaxies 
at the cores of rich galaxy clusters interact with the sur- 
rounding intracluster medium (ICM) and, in some cases, 
influence it in signiflcant and complex ways (e.g., Bohringer 
et al. 1993, 1995; Carilh, Perley & Harris 1994; Harris, Car- 
ilh & Perley 1994; Heinz, Reynolds & Begelman 1998). With 
the new generation of X-ray telescopes, Chandra and XMM- 
Newton, observations of these interactions have been more 
detailed than ever. For many systems it has been confirmed 
that there are cavities in the ICM which show up as de- 
pressions in the X-ray surface brightness and appear to be 
evacuated by the expanding radio lobes (e.g., Hydra-A; Mc- 
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Namara et al. 2000, Abell 2052; Blanton et al. 2001, Virgo- 
A; Young, Wilson & Mundell 2002, Perseus-A; Fabian et al. 
2000, 2003a). Old radio lobes which have remained intact 
through magnetic or viscous efi'ects (Reynolds et al. 2005, 
Kaiser et al. 2005, Jones & de Young 2005) are also believed 
to be responsible for the so-called "ghost" cavities seen in 
many clusters (e.g., Perseus, Fabian et al. 2003a; Abell 2597, 
McNamara et al. 2001; AbeU 4059, Heinz et al. 2002, Choi 
et al. 2004) . A second indication of the importance of radio- 
galaxy/ICM interactions is the cooling flow problem, that is 
the lack of cooled material in cluster cores despite the fact 
that the radiative cooling time of many ICM cores is rather 
short. Observations by XMM-Newton have shown that the 
gas does not in fact cool below 1-2 keV (the so-called tem- 
perature floor) begging us to identify a source of heating 
which can offset the radiative cooling. 
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At the present time, the main observational diagnostics 
for radio-galaxy /ICM interactions are the spatial mapping 
in the X-ray band of temperature, pressure, cooling time 
and metal abundances of the ICM across the interaction re- 
gion. However, with current data we lack any direct probe 
of the dynamics of the interaction, that is the velocities of 
the fluid disturbances, shocks and turbulence. A major hope 
for the Suzaku observatory was that it would open a direct 
window on the dynamics of these interactions through the 
velocity profiles of the ICM emission lines (Briiggen, Hoeft 
& Ruszkowski 2005); of course, the Suzaku cryostat failure 
and subsequent loss of the high-resolution X-ray Spectrom- 
eter (XRS) prevents this goal from being realized. Barring 
the emergence of a new mission possessing a high-resolution 
X-ray spectrometer, this important science goal must now 
await the launch of Constellation-X. 

The current baseline design for Constellation-X has 
a requirement on the spatial resolution of 15 arcsec (half- 
power diameter). Coupled with its superior spectral resolu- 
tion, this is sufficient to allow detailed mapping of radio- 
galaxy driven ICM dynamics through emission line profiles 
in the nearest systems such as the Virgo and Perseus clus- 
ters. For more distant systems, however, the interaction re- 
gion is contained within a small number of spatial reso- 
lution elements, hampering the ability of the observation 
to constrain robust dynamical signatures. In this paper, 
we discuss the complementary technique of probing these 
radio-galaxy/ICM interactions through absorption line spec- 
troscopy of the X-ray luminous core of the central radio 
galaxy, looking for jet-induced changes to the strength and 
line profile of the absorption lines present even in a static 
ICM. While detection and characterization of these absorp- 
tion line features is more observationally demanding, it has 
the major advantage of probing the kinematic state of the 
gas through a well-defined "core sample" of the galaxy clus- 
ter, even if ones X-ray observatory has only limited spatial 
resolution. We focus on predictions of resonance K-shell ab- 
sorption features of H- and He-like oxygen and iron arising 
in the ICM around an active as well as inactive radio galaxy, 
examining the velocity structure and strength of the lines. 

These radio-galaxy/ICM interactions are particularly 
interesting and important since they appear to be the lo- 
cal and accessible examples of the more general phenomena 
of "AGN feedback" and potentially crucial for some aspects 
of cosmological structure formation (e.g., see Benson et al. 
2003 and references therein). However, there is another im- 
portant reason to study this phenomenon; it is of direct rel- 
evance to our observational understanding of the Warm-Hot 
Intergalactic Medium (WHIM). The filamentary WHIM is 
believed to be the repository for half of the baryons in the 
local Universe and so is of obvious interest to cosmologists. 
It is possible to detect the WHIM by wide field imaging of 
diffuse EUV and soft X-ray emission (see for example Kaas- 
tra et al. 2003). However, the most important diagnostic of 
the WHIM are UV and X-ray absorption lines of oxygen 
seen in the spectra of bright AGN. Narrow oxygen OVI- 
OVIII absorption lines at a redshift intermediate between 
zero and that of the AGN are taken to be good candidates 
for the WHIM. However, all AGN towards which such lines 
have been claimed are radio- loud AGN (see Section 5). The 
question arises whether these absorption line features are, 
in fact, outflows associated with the jetted AGN. While the 



possibility remains that these features could be due to out- 
flows associated with the central engine itself, we demon- 
strate that it is not possible to associate WHIM-like features 
with radio-galaxy/ICM interactions. 

We predict the absorption line properties of radio- 
galaxy/ICM interactions by using the XSTAR photoioniza- 
tion code to post-process the hydrodynamics simulations of 
Reynolds, Heinz & Begelman (2002). Section 2 of this paper 
describes the underlying hydrodynamic simulations of the 
radio galaxy we use. Section 3 explains the simulation of 
our spectra using XSTAR. In Section 4, we summarize the 
properties (strengths and velocity structure) of the absorp- 
tion lines that we predict as a function of the inclination with 
which we view the radio galaxy, the age of the radio-galaxy, 
and the properties of the ICM into which the radio-galaxy 
is expanding. Section 5 argues that these interactions can- 
not produce "WHIM-like" features, thereby strengthening 
the WHIM interpretation of these intermediate redshift ab- 
sorption lines. Section 6 focuses on the ability of current 
future X-ray observatories to study these absorption lines 
and hence prove these radio-galaxy/ICM interactions. Fi- 
nally, Section 7 presents our conclusions. 



2 THE UNDERLYING HYDRODYNAMIC 
SIMULATIONS 

The hydrodynamic simulation underlying this work is the 
canonical axisymmetric simulation of Reynolds, Heinz k. 
Begelman (2002); for convenience, we briefly summarize the 
relevant details of the simulation here. 

This 2-d hydrodynamic sirrmlation explores the evolu- 
tion of an active jet in a radio galaxy as it interacts with a 
surrounding ICM, and the resulting structure after the jet 
is shut down. Initially, we start with an isothermal spherical 
gaseous atmosphere with a (number) density proflle given 
by 



The (fixed) gravitational potential is defined such that this 
atmosphere is in hydrostatic equilibrium. Back-to-back su- 
personic jets are injected in initial pressure balance with a 
density of njct = no/ 100 and Mach number (with respect to 
the internal jet sound speed) of 10. The half-opening angle 
of the jet is 15° , and the starting point for the jets is a radius 
of ro/20. As discussed in the Appendix of Reynolds, Heinz 
& Begelman (2002), this choice of parameters results in a 
cocoon in which the Kelvin-Helmholtz growth rates approx- 
imately match those calculated for real systems. Code units 
are related to physical units by parameters fixed to the back- 
ground medium. For most runs we use the following canon- 
ical conversion factors: ro = 100 kpc, cism = 1000 km s~^, 
no = 0.01 cm"^ and To = 4.4 x 10^ K. This gives a total 
kinetic luminosity of the jets of 9.3 x lO** ' erg s^^ , making 
the simulation relevant for powerful radio sources located 
in rich galaxy clusters, like Cygnus A. We will refer to this 
as our "large cluster case". Since many radio galaxies are, 
however, found in smaller and cooler galaxy clusters, we also 
consider a different scaling (ro = 10 kpc, c/sm = 500 km s~ ^ , 
no = 0.05 cm~^ and To = 1.1 x 10^ K) which results in a 
total kinetic source luminosity of 1.16 x 10**ergs~^ and is 
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Figure 1. Density map of the jet and its close environment. In 
our large-cluster ( "Cygnus-A" ) scaling , the jet in this snapshot 
has been active for 30 Myr and the region displayed is the central 
200 kpc X 300 kpc of the cluster. 



closer to objects like Hydra A and Virgo A. We will refer to 
this as our "small cluster case". In the first case, one code 
unit of time corresponds to 50 Myr. The total length of the 
simulation is 240 Myr, the jet is turned off after 50 Myr. In 
the second case, one code unit of time is 10 Myr, the total 
length is about 50 Myr, and the jet is turned off after 10 
Myr. For more detailed information on the simulations, see 
Reynolds, Heinz & Begelman (2002). As an example of the 
simulation results, a density map is given in Fig. 1. 

It is important to note that since this simulations were 
performed, Chandra observations have revealed several cases 
where the radio lobes of the central source are surrounded 
by a significant amount of cold gas. The ICM structures 
surrounding Perseus-A provide an excellent example (see 
Fabian et al. 2003a; Schmidt, Fabian & Sanders 2002). At 
the current time, the physical mechanisms underlying these 
cool rims are not understood and they cannot be reproduced 
in any first principles simulations. In systems containing 
such rims, the observed oxygen OVII and OVIII absorption 
lines towards the central AGN are likely to be dominated by 
the contribution from gas in the rim; according to Schmidt, 
Fabian & Sanders (2002), the temperature and density in 
the cold gas rim of the Perseus cluster are 3.6 X 10'^ K and 
0.12 cm~^ respectively. Using simple calculations based on 
ShuU & van Steenberg (1982), this gives an OVIII equivalent 
width of 0.15 eV, larger than that predicted from the entire 
rest of the ICM atmosphere. Therefore, our results presented 
here are valid only for systems that do not possess these cold 
rims, or for lines of sight that do not intercept such rims. 



3 SIMULATION OF THE ABSORPTION 
FEATURES 

The hydrodynamic simulation gives us the density, temper- 
ature and velocity of the ICM at every point in space as a 
function of time throughout the modelled radio-galaxy/ICM 
interaction. To simulate the absorption line spectrum that 
one would see towards the central radio galaxy, we used the 
photoionization code XSTAR (developed by Tim Kallman; 
see Kallman & Bautista 2001) to calculate the opacity of 



the ICM along radial rays that originate at the AGN. We 
assume a coUisionally-ionized plasma in thermal and ioniza- 
tion equilibrium. These assumptions are readily verified for 
the ions of interest. 

We take the temperature, density and velocity di- 
rectly from the hydrodynamic simulation. The frequency- 
dependent opacity is calculated including both thermal line 
broadening and line-of-sight Doppler shifts. We choose an 
energy resolution of AE/E — 3.29 x 10"*, giving a velocity 
resolution of about 100 km/s. 

XSTAR includes 13 elements, and the elemental abun- 
dances are important for the resulting line spectrum. We 
chose abundances measured by Allen & Fabian (1998), 
Loewenstein & Mushotzky (1996) and Baumgartner et al. 
(2005), which we summarize in Table 1. 

All publicly available photoionization codes have some 
problems dealing with very low densities and very high tem- 
peratures (> 10^ K); XSTAR is no exception. However, such 
high temperatures only occur in early times of the simula- 
tions, when the jet is still on, and then only in the inner- 
most grid cells, very close to the source. We can safely ignore 
those cells (through the inclusion of a threshold temperature 
above which the simulation cell is skipped) since very high 
temperature gas has no line absorption anyway. The only 
case where a considerable number of grid cells is affected is 
along the jet-ax;is 9 < 1° . Due to the forced symmetry of the 
simulations the region S < 1° is unphysical because material 
gets "stacked" on the jet axis. The minimum line-of-sight an- 
gle that we choose to analyze is 9 = 5.5°, appropriate for 
blazar-type sources. 



4 RESULTS 

In this section, we discuss the impact that the radio-galaxy 
activity has on the observed X-ray absorption lines in our 
simulated system. The astrophysical relevance of these re- 
sults, and the ability of current and future instruments to 
measure the predicted absorption lines will be addressed in 
the following sections. 

We performed two sets of absorption line simulations. 
The first set examines different lines of sight towards the 
AGN, at four chosen times of the jet evolution. Using the 
large cluster case, we chose to examine the spectrum of the 
system at 30 Myr (when the jet is fully active), at 50 Myr 
(just when the jet is turned off), 100 Myr (well after it is 
turned off) and at the end of the simulation at 240 Myr. Our 
second set of calculations looks at the continuous evolution 
in time, along five representative lines of sight which make 
an angle 9 — 5.5, 12, 30, 60 and 90 degrees with the jet axis. 

While our calculations actually include the thousands of 
lines that are incorporated into XSTAR, we focus our exam- 
ination on the K-shell resonant absorption lines of FeXXV 
(6.70 keV), OVII (0.574 keV) and OVIII (0.654 keV) lines 
since these are expected to be the strongest and/or most 
interesting. After examining the possibility of observing the 
features with Constellation-X, we decided to also carefully 
examine the SiXIII line (1.86 keV) since it turned out to 
be particularly strong and at an energy that is particularly 
well-tuned to Constellation-X. We note that the FeXXVI 
lines (6.95 keV and 6.97 keV) will not be discussed in much 
detail since they prove too weak to be detected even with 
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Element 


c 


N 


O 


Nc 


Mg 


Si 


S 


Ar 


Ca 


Fc 


Ni 


Z/Zq 


0.34'^ 


0.34=^ 


0.48'" 


0.62^ 


0.376'" 


0.63<: 


0.26= 


0.17<: 


0.275<: 


0.401= 


1.02= 



Table 1. Cosmic element abundances compared to solar. ''Allen & Fabian 1998, ''Loewenstein & Mushotzky 1996, ^Baumgartner et al. 
2005 



n 



energy [eV] 

Figure 2. The simulated XSTAR spectrum for 6 
t=50 Myr (large-cluster case). 



5.5° and 



our large cluster case, the FcXXV and FcXXVI equivalent 
widths are 11.65 eV and 0.086 eV (corresponding to column 
densities of A^'pexxv = 1-32 x 10^'' cm^^ and iVpeXxvi = 
2.81 X 10^^ cm"^) whereas for OVII, OVIII and SiXIII we 
predict 0.05 eV, 0.155 eV and 1.536 eV (corresponding to col- 
umn densities of Novii = 6.46 x 10^'* cm^'^, Noviii = 
1.68 X lO^'' cm-2 and Nsixiii = 1.84 x 10^*^ cm'^). 

For the small cluster case, there is no FeXXVI line 
since the ICM docs not get hot enough for appreciable 
H-like iron to exist. FeXXV and SiXIII have equivalent 
widths of 7.64eV and 11.97 eV (corresponding to column 
densities of ATpexxv = 8.67 x 10^® cm~^ and Nsixiii = 
1.43 x 10^^ cm"^), OVIII and OVII have equivalent widths of 
1.004 eV and 0.548 eV (corresponding to column densities of 
Noviii = 1-1 x lO^'' cm^^ and Novii = 7.12 x 10^" cm^^). 
The oxygen and SiXIII lines are stronger than in the large- 
cluster case since there is more cooler material to contribute 
to the absorption, but the FeXXV line is weaker due to the 
low assumed temperature of the ICM. 
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Figure 3. The strength of the most prominent lines in the sim- 
ulated spectrum for = 5.5° and t =50 Myr (large-cluster case). 

Constellation-X. For all cases, we examine the optical depth, 
equivalent width, column density and velocity shifts. Fig- 
ure 2 shows an example for the simulated full spectrum for 
9 = 5.5°, t=50 Myr (large cluster case). Figure 3 shows the 
strengths of the most prominent lines in this case. 

4.1 The undisturbed material 

Even the undisturbed ICM produces significant oxygen and 
iron line absorption (Krolik & Raymond 1988; Sarazin 1989). 
Undisturbed lines arise from material which is either ahead 
of the shock front, or which has settled down again since 
enough time has passed since the jet was turned off. For 



4.2 Disturbances introduced by the jet 

When the AGN jet moves through the ICM, it influences 
it in numerous ways. First, it can induce motions in large 
amounts of material. As the ICM gets swept up in the shock- 
front, some of it moves outwards with the shocked shell and 
some gets caught in a backflow. In the late stages of the 
evolution of the system, a significant amount of ICM falls 
back into the core regions as tfie cocoon collapses and buoy- 
antly rises into the cluster atmosphere. All of these effects 
will imprint signatures (broadening, and the development 
of blueshifted and redshifted wings and peaks) in the ab- 
sorption line profiles. The second strong effect is a change 
of the ICM temperature. Strong shocks associated with the 
early phase of jet activity will appreciably heat some regions 
of the ICM, typically raising it to temperatures sufficient to 
fully ionize oxygen. On the other hand, unshocked ICM that 
is caught in the "updraft" from a buoyant cocoon will adi- 
abatically decompress and be cooled (strengthening oxygen 
absorption features). 

Clearly, the kinematics and the thermodynamics are 
coupled. We might expect (and, as discussed below, con- 
firm) it to be appreciably harder to form high velocity oxy- 
gen absorbers than iron absorbers, since any strong shock 
responsible for accelerating the gas will inevitably heat it to 
the point where oxygen becomes fully ionized. On the other 
hand, we find iron lines which, in the strongest cases, dis- 
play a double peak structure, with one peak being at the 
rest energy of the line and one being clearly blueshifted but 
still attached to the line itself. In those cases, the blueshifted 
components arise from a large amount of shock accelerated 
material outflowing at fairly high speeds, whereas the peak 
at the rest energy comes from still undisturbed material in 
front of the jet. 
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Representative results for oxygen and iron line profiles 
and equivalent widths are shown in Figs. 4-7. It is instructive 
to separate the discussion of the cocoon inflation phase of 
activity from that of the cocoon collapse and buoyant plume 
formation phase. During the cocoon inflation phase, the ra- 
dio galaxy activity inflates an over-pressured cocoon which, 
in turn, drives a strong shock into the ICM. The high ICM 
temperatures produced by this shock leads to a rapid drop 
in the equivalent widths of the OVII, OVIII, SiXIII and 
FcXXV lines with rather little change in the line profile. The 
decrease in line equivalent widths are rather more dramatic 
for lines of sight close to the jet axis simply due to the in- 
creased path-length along which the ambient line-absorbing 
ICM has been scoured out by the cocoon shock. 

Towards the end of the cocoon-inflation phase, the co- 
coon pressure becomes comparable to that of the ambient 
ICM and the shock weakens. After that, the source transits 
into the cocoon-collapse phase in which the shocked ICM 
surrounding the sides of the cocoon (i.e., the equatorial re- 
gions with respect to the jet axis) falls back towards the 
cluster center. The result of this infall is to "squeeze" the 
cocoon plasma and transform it into two buoyantly rising, 
mushroom-shaped plumes (see Fig. 1 of Reynolds, Heinz & 
Begelman 2002). During these times, the equivalent widths 
of the OVII, OVIII, SiXIII and FeXXV absorption fines 
gradually recover to approximately their initial values. How- 
ever, the complex dynamics of the ICM during this phase 
create imprints in the line profiles. For lines of sight close to 
the jet axis, both the oxygen and iron absorption lines de- 
velop blue wings corresponding to absorption by the outward 
moving ICM shell and, latter, ICM that is being dragged out 
of the core regions of the cluster in the wake of the buoyant 
ICM plume. These features arc subtle in the case of the oxy- 
gen lines, but are rather dramatic in the CEise of the FeXXV 
lines. In all cases, the velocities characterizing the blue-wing 
are less than but of the order of the ICM sound speed. In- 
terestingly, for lines of sight that make a large angle with 
the jet-axis, the FeXXV line displays a subtle redshift of its 
centroid corresponding to the actual inward collapse of the 
ICM core (Fig. 7). 



5 IMPLICATIONS FOR THE WHIM 
OBSERVATIONS 

There has been much recent excitement about the possi- 
ble detection of the WHIM through the OVII and OVIII 
K-shell resonant absorption lines. This is of obvious impor- 
tance given that these WHIM filaments are expected to be 
the repository for half of the baryons in the local Universe. 
The most robust detection to date was obtained by a Chan- 
dra High Energy Transmission Gratings (HETG) observa- 
tion of the blazar Mrk 421 during an outburst (Nicastro et 
al. 2005a,b). Two absorber systems were detected in both 
OVII and NVII at velocities of cz = 3300 ± 300 km s"^ and 
cz — 8090 ± 300 (to be compared with the recession velocity 
of the blazar of CZ = 9000 km s'^ ). The OVII Ka EWs 
of these two systems were measured to be 0.080 ± 0.021 eV 
(3.0 ± 0.8mA) and 0.059 ± 0.021 eV (2.2 ± 0.8mA). Nicastro 
et al. (2005b) demonstrate that the column density distri- 
bution implied by these detections is consistent with the 
notion that the WHIM filaments do, indeed, balance the 



baryon budget of the local Universe. While this is the most 
robust detection of z > WHIM, it was not the first. Fang 
et al. (2002) claim a detection of the OVIII Ka line with 
an EW of 0.41 eV from a system with cz = 16, 600 km s~^ 
towards the blazar PKS 2155-304 (which has a systemic ve- 
locity of cz = 34, 800kms~^ ), although the significance of 
this detection (which is based on finding a single line in a 
blind search of the spectrum) has been questioned. McKer- 
nan et al. (2003) also claim a OVIII line with a velocity of 
cz = 4400 km s"^ and a EW of 2.17 eV towards the broad 
line radio gala^xy 3C120 (which has a systemic velocity of 
cz = 9900 kms"^ ). Other recent detections of z > WHIM 
via X-ray absorption lines have been reported towards the 
BL-Lac object PKS 0548-322 (Barcons et al. 2005) and the 
radio-loud quasar H1821-I-643 (Mathur, Weinberg & Chen 
2003). 

While observational biases are evident, it is striking that 
all of the claimed detections of z > WHIM arc along lines 
of sight towards radio-loud AGN. The question arises as 
to whether some aspect of jetted AGN can produce high- 
velocity and detached X-ray absorption lines. With our cal- 
culations, we can address one particular cleiss of models, 
i.e., those where the high- velocity X-ray absorption lines are 
produced in an jet-induced ICM outflow. In particular, we 
search for phases of activity in which one would see narrow 
and detached oxygen absorption lines. In fact, at no point 
in time do we ever see such features. Instead, any high ve- 
locity material always appears as a bluesliiftod tail on the 
main zero velocity absorption line. Furthermore, our veloc- 
ity shifts are on the order of the sound speed in our ICM. 
Any ICM that is accelerated to appreciably higher velocities 
will be shock heated to a point where it can no longer pro- 
duce oxygen absorption lines (i.e., the oxygen will be fully 
ionized). 

Thus, we conclude that the observed WHIM-like ab- 
sorption line systems cannot arise from ICM/jet interac- 
tions. This leaves two possibilities. One would be nuclear 
absorption, meaning absorption from outflowing material in 
the central nucleus of the AGN, for example from an accre- 
tion disk wind. This could be tested by looking for variability 
of the line, since the region of origin would be comparatively 
small. The other possible explanation remains the WHIM. It 
is important to obtain deep high-resolution spectra of radio- 
quiet AGN in order to assess this possibility. 



6 OBSERVING ABSORPTION LINES 
RESULTING FROM ICM-RG 
INTERACTIONS 

We have used the XSPEC spectral fitting package to in- 
vestigate the possibility of observing our simulated absorp- 
tion features with Chandra, Suzaku and Constellation- X. 
To simulate ICM absorption towards a bright radio-galaxy, 
we use our calculations to predict the absorption towards 
a power- law source with photon index F = 1.6 and flux 
3 X 10~^^ ergcm^^ s^'^. These parameters are choosen to 
mimic the radio-galaxy 3C 120 (McKernan et al. 2003). Note 
that while our liydrodynamic simulations were motivated 
by Cygnus-A, the highly absorbed nature of the nucleus of 
Cygnus-A maJtes it inappropriate as a target for such ab- 
sorption line studies. Given a predicted spectrum, we use 
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Figure 4. Line profiles and equivalent widtlis as a function of time for the small inclination case of 6 = 5.5° and the large-cluster scaling. 
Upper left panel: Energy profiles of the OVIII Ka resonance absorption line. Upper right panel : Total equivalent width of the OVIII 
Ka resonance absorption line as a function of time. Lower left panel: Energy profile of the FeXXV Ka resonance absorption line. Lower 
right panel: Total equivalent width of the FeXXV absorption lines as a function of time. Note the growth of an extended blue-shifted 
wing at the end of the cocoon-inflation phase. 



the "fakeit" tool in XSPEC to simulate spectra using the 
appropriate response matrices for the high-energy transmis- 
sion gratings (HETG) on Chandra, the now in-operable X- 
ray spectrometer (XRS) on Suzaku and the microcalorimeter 
on Constellation-X. We chose the optimistic case of integra- 
tion times of 200 ks. 

To determine the detectabihty of the hnes of interest, 



we first fitted the simulated data with a powerlaw. We then 
search for the absorption lines by adding a simple negative 
Gaussian model and applying the F-test to the resulting im- 
provement in the goodness of fit. This procedure reveals that 

neither C/iandra/HETG nor Suzaku/yiYiS are able to detect 
the iron and oxygen absorption features for lines-of-sight 
close to the jet axis, because along those lines-of-sight the 
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Figure 5. Line profiles and equivalent widths as a function of viewing angle at the time that the jet shuts off (t = 50Myr) for the 
large-cluster scaling. Upper left panel: Energy profiles of the OVIII Ka resonance absorption line. Upper right panel : Total equivalent 
width of the OVIII Ka resonance absorption line as a function of inclination. Lower left panel: Energy profile of the FeXXV Ka resonance 
absorption line. Lower right panel: Total equivalent width of the FeXXV absorption lines as a function of inclination. Note the growth 
of an extended blue-shifted wing at small inclinations. 



predicted features are also the weakest. The only detectable 
feature at 6 = 5.5° is the SiXIII line. The simulated line 
as detected by Suzaku is shown in Fig. 8. Detections of the 
other features with Suzaku are possible at higher inclina- 
tions. As discussed in Section 4 however, the ICM dynam- 
ics imprint only subtle features on the line profile for these 
higher inclinations. Thus, we conclude that studies of ICiM 



dynamics using these absorption lines are not possible with 
the current generation of instruments. 

On the other hand, NASA's proposed Constellation- 
X observatory will have such improved sensitivity that 
it will easily allow the study of jet-driven ICM dynam- 
ics using absorption line profiles. Looking at the case of 
9 = 5.5° and t=50 Myr, the most prominent lines picked 
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Figure 6. Line profiles and equivalent widths as a function of time for the small inclination case oi 9 = 5.5° and the small-cluster scaling. 
Upper left panel: Energy profiles of the OVIII Ka resonance absorption line. Upper right panel : Total equivalent width of the OVIII 
Ka resonance absorption line as a function of time. Lower left panel: Energy profile of the FeXXV Ka resonance absorption line. Lower 
right panel: Total equivalent width of the FeXXV absorption lines as a function of time. Note the growth of an extended blue-shifted 
wing at the end of the cocoon-inflation phase. 



up by Constellation-X are SiXIII, FeXXV, SXV, MgXI, 
FeXXIV, NelX, OVIII, NeX and OVII. As shown in Fig. 9, 
FeXXV displays the double peaked structure with one 
strongly blueshifted second peak introduced by the out- 
flowing shocked ICM shell. SiXIII and SXV both display 
blueshift wings making a Gaussian not an adequate model 
anymore. All lines display blueshifted peaks with At; = 10® 



to 10*^ cm s ^ . The simulated FeXXV line detection is shown 
in Fig. 9. 



7 CONCLUSIONS 

In this paper we present spectral simulations of ICiM absorp- 
tion features arising from radio-galaxy /ICM interactions. 
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Figure 7. Line profiles and equivalent widths as a function of time at a viewing angle of 60° for the large-cluster scaling. Left panel: 
Energy profiles of the OVIII Ka resonance absorption line. Right panel: Energy profile of the FeXXV Ka resonance absorption line. Note 
the subtle redshift of the centroid of both lines associated with the lateral collapse of the cocoon, at t=36 Myr (solid line) . 
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Figure 8. Simulated Suzaku detection of the SiXIII line in the 
large-cluster case for an inclination oi d = 5.5°, and a time of 
t = 50 Myr. See text for details of the simulated source and ob- 
servation. 
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Figure 9. Simulated Constellation- X detection of FcXXV in the 
large-cluster case for an inclination oi = 5.5°, a time of t = 
50 Myr. See text for other details of the simulated source. Note 
the double-peak structure. The blueshifted peak corresponds to 
material in the approaching shocked-shell of ICM. 



These features can be used as a new tool to probe inter- 
actions between the ICiM and radio galaxies. The energy in- 
troduced by the jet weakens the lines already present duo to 
the static ISM and, more importantly, introduces a velocity 
structure in the form of wings and blueshifted line peaks in 
the line profile. We specifically looked at OVIII, OVII and 
FeXXV. The results were quantified as equivalent widths, 
column densities and velocity shifts of those absorption lines. 
We used the results for the oxygen lines to test the possi- 
bility of jet-induced ICM outflow giving rise to absorption 
lines which could be mistaken as WHIM in the claimed ob- 



servations. Our results largely rule that out and strengthen 
the WHIM interpretation. We also tested the possibility of 
actually observing those features introduced by the jet with 
the HETG on Chandra, the now inoperable XRS on Suzaku 
and the microcalorimeter on Constellation-X. Chandra is not 
sensitive enough to detect the features. Suzaku would have 
been barely sensitive enough to detect SiXIII and FeXXV, 
but only along lines of sight fairly far away from the jet axis, 
where the jet influence is minimal and therefore it is not pos- 
sible to learn much about the interaction dynamics between 
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the AGN and the ICM. Constellation- X will have a vastly 
superior sensitivity and is therefore much more suitable for 
that task. In our simulations, Constellation-X readily de- 
tects these features, can detect double-peaked FeXXV lines 
resulting from the expanding shocked ICM shell, and can 
detect the re-collapse of the inner ICM after the AGN hcis 
turned off. We suggest that such observations may be a pow- 
erful new diagnostic of the dynamics of radio-galaxy /ICM 
interactions. 
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